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ABSTRACT 

Motivated by the recent discovery of 30 new millisecond pulsars in Terzan 5, made using the Green 
Bank Telescope's S-band receiver and the Pulsar Spigot spectrometer, we have set out to use the 
same observing system in a systematic search for pulsars in other globular clusters. Here we report 
on the discovery of five new pulsars in NGC 6440 and three in NGC 6441; each cluster previously 
had one known pulsar. Using the most recent distance estimates to these clusters, we conclude that 
there are as many potentially observable pulsars in NGC 6440 and NGC 6441 as in Terzan 5. We 
present timing solutions for all of the pulsars in these globular clusters. Four of the new discoveries 
are in binary systems; one of them, PSR J1748-2021B (NGC 6440B), has a wide (Pb = 20.5 d) 
and eccentric (e = 0.57) orbit. This allowed a measurement of its rate of advance of periastron, 
(ju = 0.00391(18)°yr~'^. If due to the effects of general relativity, the total mass of this binary system 
is 2.92 ± 0.20 Af0 (1 a), implying a median pulsar mass of 2.74 ± 0.21 Mq. There is a 1% probability 
that the inclination is low enough that pulsar mass is below IMq, and 0.10% probability that it is 
between 1.20 and 1.44 Mq. If confirmed, this anomalously large mass would strongly constrain the 
equation of state for dense matter. The other highly eccentric binary, PSR J1750— 37A, has e — 0.71, 
and io = 0.0055(3)°yr~^, implying a total system mass of 1.97 ± 0.15 M© and, along with the mass 
function, maximum and median pulsar masses of 1.65 and 1.26 Mq respectively. 
Subject headings: binaries: general — stars: neutron — pulsars: general — globular clusters: individ- 
ual (NGC 6440, NGC 6441) — equation of state 
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1. INTRODUCTION 

A recent survey of the globular cluster (GC) Terzan 5 
with the S-band (1.7 to 2.6 GHz) receiver of the Green 
Ba nk Telescope (GBT) and the Pulsar Spigot spectrome- 
ter ([Kaplan et al. 20051 ) has met with extraordinary suc- 
cess, finding 30 new millisecond pulsars (MSPs). The 
first 21 discoveries were announced in Ransom et al. 
(2005). More recently, nine more pulsars were discov- 
ered^, including PSR J1748-2446ad, th e fastest spin- 
ning pulsar known ( Hessels et al. 200ti). These dis- 
coveries show that this observing system (henceforth 
GBT/S/PS) has an unprecedented sensitivity to MSPs 
outside the Arecibo sky, due in great part to the fac- 
tor of 3 higher gain of the GBT as compared to the 
Parkes telescope. Perhaps even more importantly, the 
higher observing frequency (~2GHz), large observing 
bandwidth (600 MHz), and relatively fine frequency res- 
olution (~0.78MHz) retain high sensitivity and time res- 
olution for MSPs with large dispersion measures (DM). 
The discovery of the 1.39 ms PSR J1748-2446ad at a 
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DM of 235.6 cm^'^pc dramatically demonstrates the sen- 
sitivity of the GBT/S/PS system to fast-spinning pulsars 
at high DMs. This motivated us to search for MSPs in 
other promising GCs with the GBT/S/PS system. 

1.1. Cluster selection 

In the high stellar density environments in the cores 
of GCs, binary systems with at least one main sequence 
(MS) star are occasionally disrupted by the intrusion of 
a neutron star (NS). The most likely outcome is an "ex- 
change encounter" , where the lightest component of the 
previous binary is ejected and a NS-MS binary remains. 
The resulting NS-MS systems can become much more nu- 
merous than those formed from evolution of primordial 
binaries. With the evolution of the MS star, it can hap- 
pen that matter will accrete onto the NS during a Low- 
Mass X-ray Binary (LMXB) phase, particularly among 
the more compact systems. When the LMXB accretion 
stops, we are left with an MSP orbiting a white dwarf 
(WD) star. 

Pooley et al. (2003) showed a strong correlation be- 
tween the number of X-ray point sources in a GC (which 
are generally compact binaries) and the rate of stellar en- 
counters within the half-mass radius, a parameter that 
they designate F. This correlation is stronger than for 
any other individual GC parameter, and applies to ac- 
tive LMXBs as well (although with only a dozen systems 
in the Galactic GC system, the statistics are poorer). P 
might therefore be a good indicator of the present rate 
of MSP formation. Pooley et al. (2003) provide a simple 
way of calculating the stellar interaction rate at the core 
of any GC: 

Tccxpj-V^, (1) 

where po is the central mass density and Tc is the core 
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TABLE 1 
The six globular clusters with highest Fc 



Globular 


£,a 


PO-^ 


Tc 


DM'^ 


Tscatt,2GHz'* 


A^psr" 


ATf 


DM Ref. 


Cluster 


(kpc) 


(log [L0 pc-3]) 


(%) 


(cm~^pc) 


(MS) 




(hr) 


g 


NGC 6388 


10.0 


5.34 


9.5 


318 


1.4 





2.2 




Terzan 5 


10.3 


5.06 


8.6 


(234.3-245.6) 


6 


3 +30 


8.0 


(1) 


NGC 6441 


11.7 


5.25 


8.0 


(230.1-234.4) 


2 


1 +3 


5.3 


(2) 


NGC 6440 


8.4 


5.28 


6.4 


(219.4-227.0) 


5 


1 +5 


8.4 


(3) 


NGC 6266 


6.9 


5.14 


5.1 


(113.4-115.0) 


4 


6 


7.0 


(4) 


Liller 1 


9.6 


5.53 


4.2 


771 


~ 3 X 102 





6.2 





'^ Distance from the Solar System from Harris (1996). These values have signifi- 

cant systematic uncertainties.*^ Central density of GC; see last column of third table in 
http://physwww.mcmaster .ca/$\sim$harris/mwgc.dat '^ For GCs with known pulsars, we present 
the DM range tor the known pulsars, with references in the last column, otherwise we predict the DM for 
the GCs estimated distance using the Cordes & Lazio (2001) model. These can be off by a factor of up 
to about 2.^* Extrapolated from other frequencies using Tscatt,2GHz = ''"scatt(^)(i^/2GHz)^'*. In the case 
of Terzan 5, the scattering was measured by Nice & Thorsett (1992) as 700 /xs at 685 MHz. Otherwise, 
we scale the scattering time at IGHz predicted by the Cordes & Lazio (2001) model for the DM in the 
previous column. These can be off by up to an order of magnitude."^ Number of previously known pulsars 
+ number of new pulsars.' Amount of time the GC is visible at the GBT per day^ (1) Ransom et al. 
2005, (2) Possenti et al. 2006, this work, (3) Lyne, Manchester & D'Amico 1996, this work, (4) Possenti 
et al. 2003, Chandler 2003. 



TABLE 2 
Parameters for three globular clusters 



NGC 6388 



NGC 6440 NGC 6441 



Right Ascension of center, a . . . . 17*^36™ 17;0 

Declination of center, <5 -44°44'06" 

Galactic longitude, I (°) 345.56 

Galactic latitude, b (°) -6.74 

Cluster distance, D (kpc) 10.0 

Angular core radius, dc (arcmin) 0.12 

Half-mass radius, dy^ (arcmin) . . 0.67 

Metallicity (log[Fe/H]) -0.60 

v^{0) (km s-l) 19.13*^ 

ac^ (10-9ms-2) -0.46 



7h48™52?7 


17^50^12=9 


20°21'37" 


-37°03'05" 


7.73 


353.53 


3.80 


-5.01 


8.2'' 


IS.S'' 


0.13 


0.11 


0.58 


0.64 


-0.34 


-0.53 


13.01'' 


i9.5i^:t^ 


-0.029 


-0.46 



Note. — Parameters as in Harris (1996), except where indicated. 
'^ Valenti et al. (2007).'' Webbink (1985). '^ Dubath et al. (1997).'' Calculated 
using the data above and the Kuijken & Gilmore (1989) model of the Galaxy. 



radius. This parameter is perhaps more relevant than F 
for the formation of LMXBs because mass segregation 
constrains neutron stars to within a few core radii of the 
center of the GC. 

We calculated Tc for all Galactic GCs with the cen- 
tral luminosity densities and core radii given in Har- 
ris (1996)^, and assuming a constant conversion factor 
between central luminosity density and central mass den- 
sity. The six best GCs are listed in Table[TJ with their Tc 
presented as a percentage of the total T^ summed over 
all Galactic GCs. This table represents our initial sur- 
vey expectations; we hypothesized that the percentages 
should be closely related to the fraction of the total GC 
MSP population contained in each GC. 

For the first phase of the survey, we selected the GCs 
in Tabic [1] that had not been previously searched with 
the GBT/S/PS system: NGC 6388, NGC 6440 and 
NGC 6441. Liller 1 was searched at the same time as 
Terzan 5 by Ransom et al. (2005) and NGC 6626 has 
been searched by Jacoby et al. (2002). Despite the large 

^ See |http:/ /physwww.mcmaster.ca/~harris/mwgc.dat for an 
updated TiSTofTJU parameters. ' 



predicted stellar interaction rates, only two pulsars were 
known in these GCs: PSR B1745-20A (NGC 6440A; 
Lyne et al. 1996) and PSR J1750-37A (NGC 6441A; 
Possenti et al. 2006). The spin periods and DMs of 
NGC 6440A and NGC 6441A are 288.6ms and 111.6ms 
and 220 and 234 cm~'^ pc, respectively. We assumed that 
the lack of other known pulsars in these high-DM GCs 
was largely due to a bias of the earlier surveys against the 
detection of faster-spinning MSPs. The high central fre- 
quency and improved time resolution of the GBT/S/PS 
system greatly diminish pulse broadening due to multi- 
path propagation and dispersive smearing, and would 
therefore alleviate the bias against the detection of other 
MSPs in these GCs. 

1.2. Data taking 

From 2005 April 18 to August 21, we observed 
NGC 6388, NGC 6440 and NGC 6441 on three different 
occasions each, using the GBT's S-band receiver with 
the Pulsar Spigot. The longest integration times used 
for each GC are very close to their maximum track times 
indicated in Table [51 

These observations take advantage of the relatively 
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radio-frequency-interference (RFI) clean band between 
1650 and 2250 MHz. The Pulsar Spigot is built around 
the GET correlati on spectrometer an d is available in a 
number of modes (jKaplan et al. 20051 ). For most obser- 
vations, we used the Spigot in mode 2 where two orthog- 
onal polarizations with a total bandwidth of 800 MHz 
are digitized with three levels, autocorrelated with 1024 
lags, and integrated for 81.92 /xs. The resulting autocor- 
relation functions are summed and then written to disk 
as 16-bit integers. After Fourier transforming the lags, 
we have 768 usable channels across the 600-MHz band- 
pass (0.78125-MHz channels). At a frequency of 2 GHz 
and a DM of 220 cm~^ pc, the dispersive smearing across 
one channel is 178 /iS. Recently, we have been using mode 
14, this mode writes the 81.92-/is integrations to disk as 
8-bit integers, but with twice the frequency resolution, 
resulting in a dispersive smearing across one channel for 
a DM of 220 cm'^ pc of 89 /zs. 

In the case of NGC 6440 and NGC 6441, the discovery 
of new pulsars lead us to observe these two GCs on a 
total of 37 and 24 extra occasions respectively, from 2005 
November 5 to 2007 March 28 (Spigot mode 2), and then 
from 2007 October 10 to November 1 (Spigot mode 14), 
the latter data has only been added to the ephemeris 
of NGC 6440B (see fj3]). The typical observing session 
spends the first 2-3 hours on NGC 6440, then observes 
NGC 6441 for 2.5 to 4 hours, and then spends 2 to 3 
hours again on NGC 6440. The total usable observing 
times are about 90-100 hr for each pulsar in NGC 6440 
and NGC 6441. 

On two occasions, we observed NGC 6440 and 
NGC 6441 at 820 MHz with a bandwidth of 50 MHz in 
order to determine the DMs of the pulsars more accu- 
rately and to search for pulsars with spectral indices too 
steep to be detected at 1950 MHz, or that lie outside the 
area covered by the S-band beam. Because the telescope 
beam widths at these frequencies (6f5 at 1.95 GHz and 
15' at 820 MHz) are significantly larger than the half- 
mass radii of these GCs (see Table [2|) , we observed only 
single positions centered on the GCs. 

2. SEARCH AND DISCOVERY OF NEW PULSARS 

After masking the interference (in the manner de- 
scribed in detail in Hessels et al. 2007), the search data 
(taken in Spigot mode 2) were dedispersed in 48 sub- 
bands over the whole 600-MHz usable band at the "av- 
erage" DMs for NGC 6440 (223 cm'^ pc) and NGC 6441 
(233 cm~^ pc), reducing the volume of data to be stored 
by a factor of ~20. For each sub-band set we created 
~40 dedispersed time series covering a range of trial DMs 
iblOcm^'^pc around the nominal GC DMs in steps of 
0.5 cm~^ pc. If a pulsar has a DM halfway between two 
trial DMs, that introduces an extra smearing of 168/iS, 
resulting in a total maximum smearing of 245 fis. We 
also barycentered each time series using TEMPO^". 

We Fourier transformed each time series and searched 
them using a frequency domain acceleration search tech- 
nique (Ransom et al. 2002) that allows a fully coherent 
search in frequency/frequency-derivative space as well as 
the incoherent addition of a number of harmonics (1, 2, 
4, 8 or 16). The technique has an optimal sensitivity for 
pulsars in binary systems with an orbital period much 

^^ http://www.atnf.csiro.au/rcscarch/pulsar/tcnipo/ 



longer than the observation time (Porb > lOTobs)- For 
this reason, the acceleration searches were performed on 
all the observations for the entire length of each observa- 
tion as well as on overlapping segments of various shorter 
durations (typically 27 and 82 minutes). A more detailed 
description of the search procedure is presented in Hes- 
sels et al. (2007). 

We discovered eight new MSPs, all in S-band data, five 
of them in NGC 6440 and three in NGC 6441. Prehm- 
inary results on these objects were presented in Begin 
(2006). Their pulse profiles, together with those of the 
previously known objects, are presented in Figure [TJ 
These pulse profiles were obtained by adding all detec- 
tions of the pulsars made in mode 2, carefully excluding 
all sub-integrations in time or frequency that are cor- 
rupted by significant RFI. 

We found no pulsars in NGC 6388, possibly due to 
dispersive smearing and/or interstellar scattering being 
substantially higher than the prediction given in Table [TJ 
The search for pulsars in NGC 6388 was further hin- 
dered by the small time the GBT can observe the GC 
(2.2 hrs; see Table [T]). Furthermore, many of the pulsars 
in Terzan 5, NGC 6440 and NGC 6441 were found by 
searching follow-up timing observations, none of which 
were made for NGC 6388. 

2.1. Flux densities 

The 1950 MHz flux densities (S'1950) of all the pulsars in 
NGC 6440 and NGC 6441 and their pseudo-luminosities 
at the same frequency (^1950 = »5'i95o-D'^) are listed in 
Tables[3] (isolated pulsars) and|4] (binary pulsars). These 
were determined from the pulse profiles in Fig. [1] by as- 
suming that the off-pulse r.m.s. is as predicted by the 
radiometer equation for the total observing time. In cal- 
culating the system's equivalent flux density, we used a 
system temperature of 26K towards NGC 6440 and 28 K 
towards NGC 6441, a gain of 1.95 K/Jy for the GBT at 
S-band and multiplied the result by 1.235 for the 3-level 
efficiency factor of the correlator. 

These flux densities indicate that some of the new dis- 
coveries, particularly in NGC 6441, are as faint as the 
faintest pulsars in Terzan 5, suggesting that we have 
reached a similar sensitivity in our survey as in that re- 
ported by Ransom et al. (2005). 

We can compare the sizes of the pulsar populations in 
these clusters with that of Terzan 5 by comparing the 
distribution of their pseudo- luminosities. The observing 
system and observing frequency are the same, the DMs 
of the pulsars are very similar, the scattering times are 
similarly small, the length of the integrations are similar 
(Table[T]) and the search software is the same, making the 
search similarly sensitive to fast-spinning pulsars (i i2.2p . 
Furthermore, the software used to make the flux density 
measurements is similar (but not exactly - some refine- 
ments have recently been added to the program). All 
these factors make the flux densities directly comparable; 
the pseudo-luminosities will then be directly comparable 
as well if the GC distances are accurate. 

We have fitted a luminosity function of the type A^(> 
L) = kiL~^-^ (derived by Hessels et al. 2007 as a 
best fit to the luminosities of 41 isolated pulsars in M5, 
M13, M15, M28, Terzan 5, 47 Tucanae, NGC 6440 and 
NGC 6441 above L1400 = 1.5 mJy kpc^), where iV(> L) 
is the number of pulsars above a pseudo-luminosity of L 
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Fig. 1. — Average 1950 MHz pulse profiles for the 10 pulsars known in the GCs NGC 6440 and NGC 6441, obtained by adding all 
observations made with mode 2 of the Pulsar Spigot. These profiles cover one full rotation. The pulsars in binary systems have an asterisk 
after their names, and the newly discovered pulsars have their names in boldface. The horizontal width of the rectangles indicates the 
system's total time resolution in mode 2, including the effects of dispersive smearing, relative to each pulsar's spin period. The dip in power 
after the main pulse of NGC 6440A is likely a Spigot artifact. 
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Fig. 2. — Cumulative distribution of pseudo-luminosities at 1950 
MHz for the pulsars in NGC 6440, NGC 6441 and Terzan 5. The 
pseudo-luminosities are estimated assuming the most recent GC 
distance estimates (Ortolani et al. 2007, Valenti et al. 2007). 



and ki is the number of pulsars brighter than 1 niJy kpc^. 
If we assume the distances Usted in Harris (1996, see Ta- 
ble [J), we obtain ki = 32.8,7.9 and 5.7 for Terzan 5, 
NGC 6440 and NGC 6441 respectively. This would indi- 
cate that Fc does a poor job at predicting the number of 
MSPs in these GCs. 

However, recent studies of these GCs have changed 
these distance estimates (see Table [2]) , sometimes dras- 
tically: in the case o f Terzan 5, the most recent distance 
estimate is 5.5 kpc (jOrtolani et al. 2007T ). It should be 
noted that these new distance estimates still have large 
systematic uncertainties (Heinke and Ortolani, private 
communication). These values affect our prediction for 
Tc- the total absolute visible luminosity of the core of 
the GC is proportional to D^ , but the volume of the core 
is proportional to D^, therefore po oc D^/D^ = D~^. 
Inserting this in Eq. [l] and keeping in mind that Vc is 



also proportional to D, we obtain Tc oc D^^^. With the 
revised distance estimates, we obtain revised Tc values of 
6.3, 6.3 and 8.6% for Terzan 5, NGC 6440 and NGC 6441 
respectively. 

When we use the new distance estimates to re-calculate 
the pseudo-luminosities of the pulsars in these three GCs 
(Fig. [2]), we obtain ki = 10.6,7.5 and 7.3 for Terzan 5, 
NGC 6440 and NGC 6441 respectively. Given the uncer- 
tainties in the distance estimates, this is consistent with 
all the clusters having similarly sized pulsar populations, 
and Fc being a good predictor of the number of pulsars 
in a GC. 

If the fci numbers reflect the sizes of the pulsar popula- 
tions, then assuming the Hessels et al. (2007) luminosity 
law, we should detect 
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(2) 
pulsars in NGC 6440 and 4.1 pulsars for NGC 6441". 
For NGC 6441, we detect as many pulsars as expected 
from the luminosity distribution given its greater dis- 
tance, but we don't detect faint (10 - 17 /iJy) pulsars 
in NGC 6440. This could happen if some of them arc in 
tight binary systems, which would make them more diffi- 
cult to detect. In any case, the most likely reason why we 
detect more pulsars in Terzan 5 is simply because that 
GC is closer to us. 

To summarize, we can see from the discussion above 
that the distance to the GC does not have a major im- 
pact on the estimate of Tc (Tc oc D^^^), but it has a fun- 
damental impact in estimating pulsar populations from 
the observed flux densities (fci oc _D^^°-^). For Terzan 5, 
NGC 6440 and NGC 6441, the use of the latest distance 
estimates makes Tc a good predictor of the relative sizes 
of the pulsar populations. 

2.2. Spin period distribution 

Looking at the spin periods in Table [51 and comparing 
them with the spin periods of the pulsars in other GCs^^, 

^^ The number of pulsars in Terzan 5 we use in this calculation 
is 30. This many were detected using the same search methods 
employe d in the present survey. 

^^ See|http://www.naic.edu/'^pfreire/GCpsr.html| 
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we notice a conspicuous absence of fast-spinning objects. 
Of the 33 known pulsars in Terzan 5 (a GC for which the 
discovery observations have very nearly the same time 
resolution as in our survey), 15 have spin periods shorter 
than that of NGC 6440F, the fastest MSP we found, and 
23 have spin periods shorter than that of NGC 6441D. 

This could be a result of interstellar scattering being 
worse for NGC 6440 and NGC 6441 than it is for Terzan 
5. Looking at Tabled we can see that both NGC 6440 
and NGC 6441 have predicted scattering times lower 
than those of Terzan 5. An inspection of the pulse pro- 
files in FigurelU particularly of the shorter-period pulsars 
such as NGC 6440C and NGC 6441D, does not reveal the 
typical signature of interstellar scattering, an exponen- 
tial decay of intensity after the main peak, as one would 
expect from the small scattering timescales in Table [TJ 
The difference, if statistically significant, is likely to be 
real, not due to scattering or dispersive smearing. 

A two-sided Kolmogorov-Smirnov (2K-S) test 
([Press et al. 19921 ) gives an 87% probability that the 
spin period distributions NGC 6440 and NGC 6441 are 
those of a common population. Comparing these two 
with Terzan 5, we obtain, respectively, 1.1% and 1.8% 
probabilities for the same hypothesis. Comparing these 
two with the pulsar population in 47 Tucanae, we obtain 
probabilities of 0.13% and 0.27% respectively. 

This phenomenon is not restricted to NGC 6440 and 
NGC 6441. A recent 1.4 GHz GC pulsar survey made 
with the Arecibo telescope (Hessels et al. 2007) has 
searched for pulsars in M15. Despite being even more 
sensitive to short-period pulsars than the Terzan 5 sur- 
vey, it found none. The spin period distribution of the 
M15 pulsars detected by Hessels (2007) is indistinguish- 
able from those of NGC 6440 and NGC 6441 (2K-S prob- 
abilities of 62% and 65%), but different from that of 
Terzan 5 (2K-S probabihty of 2.6%). 

2.3. Dispersion measure distribution 

The DMs of the known pulsars in NGC 6440 range 
from 219.4 to 226.9 cm^'^ pc, a span of 7.5cm"'^pc. 
In Terzan 5, the known pulsars have DMs be- 
tween 234.3 to 243 .5 cm~'^pc, a span of 9.2cm~^pc 
(jRansom et al. 20051 ). However, two of the pulsars in 
Terzan 5 (A and D) have anomalously high DMs. Each 
of these pulsars is spatially more distant from the center 
of the GC than the other pulsars. All the other pul- 
sars have DMs between 234.32 and 239.82 cm-^pc. If 
we randomly choose six pulsars in Terzan 5, there is a 
66% probability that the set will include none of the two 
outliers. In that case, the maximum DM span will be 
5.5cm~^pc, lower than but comparable to what we ob- 
serve in NGC 6440. 

This is consistent with the idea that the variations 
in DM between pulsars are generally due to differences 
in the Galactic electron column density ^^. Since both 
NGC 6440 and Terzan 5 are at similar DMs, and have 
their pulsars spread over roughly similar areas of the sky, 
the differences in DM between the pulsars should be sim- 
ilar. The spread in DMs seems to be linearly correl ated 
with the absolute DM of the GC (jFreire et al. 2005f) . an 



^^ The exception is 47 Tucanae lIFreire et al. 200ll ). wliere the 
DM differences between pulsars are predominantly due to gas in- 
trinsic to the GC. 



idea which is supported by the observed spread of DMs 
in NGC 6440. 

For NGC 6441, where only four pulsars are known, 
the DMs range between 230.09 and 234.39cm-3pc. The 
slightly more compact distribution of the smaller number 
of pulsars is likely to produce a smaller range of DMs 
than observed in NGC 6440. 

3. TIMING 

Diffractive scintillation causes a modulation of the in- 
tensity of the pulsed signal in time and frequency. "Scin- 
tles" are the regions of the dynamic spectrum where the 
signal is amplified relative to average. The high DMs of 
these pulsars imply that the scintles have narrow band- 
widths (of the order of -20 kHz and -60 kHz at 1.95 
GHz for NGC 6440 and NGC 6441, Cordes and Lazio 
2001), and are therefore averaged out over our wide ob- 
serving band. This means that the flux densities of 
the pulsars are very steady, resulting in the very high 
detection rate during our timing campaign: even the 
faintest new MSPs are always detectable. These detec- 
tions were cross-correlated with a synthetic profile, ob- 
tained from fitting a minimal set of Gaussian curves to 
the summed profiles presented in Figure [TJ This cross- 
correlation is done in the Fourier domain, as described 
in Taylor (1992), and from it we derive pulse times of 
arrival (TO As). For the faint, isolated pulsars we calcu- 
lated about 1 TOA per observation, but for the binary 
systems, like NGC 6440B and particularly NGC 6440D, 
a significantly larger number were calculated in order to 
retain orbital phase information. 

These TOAs were fitted to a pulsar model contain- 
ing spin (its frequency v and its first derivative v) and 
astrometric (right ascension a and declination S) pa- 
rameters, plus the DM using TEMPO. The parame- 
ters obtained from these fits are presented in Tables [3l 
and [H To account for the motion of the telescope rela- 
tive to the barycenter of the sola r system, we u sed the 
DE405 Solar System ephemeris ([Standish 1998f) . The 
DMs are calculated by comparing pulse TOAs at 820 
and 1950 MHz, but they assume that there is no sig- 
nificant longitudinal evolution of the pulse profile with 
frequency. For the binary pulsars, the orbital param- 
eters (orbital period P^; semi- major axis of the pul- 
sar's orbit a projected along the line-of-sight in light- 
seconds, X = asini/c; eccentricity e; longitude of 
periastron cu; time of passage through periastron Tq; 
and in two cases the rate of advance of periastron w) 
were also fitted (see Table [^ using the the Damour 
fc Deruelle orbital mod el ([Damour fc Deruelle 1985[ : 
iDamour fc Deruelle 198^ ). The reference epoch for all 
solutions is MJD = 54000 (2006 September 9). The tim- 
ing solution of NGC 6440B is the only one that includes 
the "mode 14" data taken in October 2007; all the others 
include only mode 2 data, which ends in 2007 March 28. 

The 3-level sampling of the Spigot introduces system- 
atics on the profiles that are not well understood. Fol- 
lowing standard pulsar timing practice in such cases, we 
compensate for this by increasing the uncertainties on 
the TOAs by the small factor indicated in Tables [31 and 
[Has "EFAC", such that the reduced x^ = 1. In all cases 
except that of NGC 6440C, these factors are small, in- 
dicating that the timing models in Tables [31 and [H do a 
good job of predicting the rotational phases of the pul- 
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TABLE 3 
Parameters for the isolated pulsars in NGC 6440 and NGC 6441 



PSR B1745-20A PSR J1748-2021C PSR J1748~2021E PSR J1750-3703C PSR J1750-3703D 

(NGC 6440A) (NGC 6440C) (NGC 6440E) (NGC 6441C) (NGC 6441D) 

Observation and data reduction parameters 

^1950'' (mJy) (137 0044 (1023 OUTB OUTO 

ii95o'' (mjy kpc2) 25 3.0 1.5 2.7 1.8 

Number of TOAs 107 161 61 82 68 

Residual rms (^s) 452 38 21 100 53 

EFAC 1.72 3.37 1.02 1.17 1.25 

Reference Epoch (MJD) 54000 54000 54000 54000 54000 

Timing parameters 

a (J2000) 17'^48™52!689(2) 17'M8™51! 17320(15) 17'M8™52!80040(14) 17''50™13! 4541(7) 17''50™13!0969(4) 

<5 (J2000) -20°21'39('7(5) -20°21'53'.'81(4) -20°21'29'.'38(3) -37°03'05'.'58(3) -37° 03'06y 368(17) 

u (Hz) 3.464969947246(19) 160.59270990836(6) 61.48547652886(2) 37.63830424257(6) 194.55492232569(15) 

i> (10-15 Hz s~i) -4.7944(12) 1.543(4) -1.1812(15) 1.411(4) -18.652(9) 

DM(cm-^pc) 219.4(2) 226.95(6) 224.10(4) 230.67(2) 230.09(17) 

^ The uncertainties in these fluxes are of the order of 10%.'' In this and the next table, the distances used to calculate this parameter are 
those given in Table [2] 



TABLE 4 
Parameters for the binary pulsars in NGC 6440 and NGC 6441 



PSR J1748-2021B 
(NGC 6440B) 

Observation and data reduction parameters 

5i950 (mJy) 0047 

Ligso (mJy kpc^) 3.2 

Number of TOAs 1051 

Residual rms {/is) 46 

EFAC 1.20 

Reference Epoch (MJD) 54000 

Timing parameters 

a (J2000) 17''48™52!95290(8) 

(5 (J2000) -20°21'38'.'86(19) 

u (Hz) 59.665418222544(17) 

!> (10-15 Hz s~i) 1.1717(6) 

DM(cm-3pc) 220.922(11) 

Pb {d) 20.5500072(6) 

X (1-s) 4.466994(6) 

To (MJD) 54005.480292(7) 

e 0.5701606(15) 

a; (°) 314.31935(13) 

dj (° yr-i) 0.00391(18) 

Derived parameters 

/(M0) 0.0002266235(9) 

M(Mq) 2.92(20) 

Mp,max(A4'Q) 3.25 

M,,^in(Mo)=' 0.11 



PSR J1748-2021D 
(NGC 6440D) 



PSR J1748-2021F 
(NGC 6440F) 



PSR J1750-37A 

(NGC 6441A) 



PSR J1750-3703B 
(NGC 6441B) 



0.075 
5.0 
946 
42 
1.33 

54000 



0.017 

1.1 

52 

36 

1.35 

54000 



0.059 
10.8 
129 
86 
1.08 

54000 



0.037 
6.7 
277 
116 
1.55 

54000 



17''48'"51!64665(7) 

-20°21'07'.'414(18) 

74.097014488305(13) 

-3.2216(9) 

224.98(3) 

0.2860686769(4) 

0.397203(3) 

54000.1053987(3) 

[0.0] 





17'' 48™52! 3339(3) 

-20°21'39'.'33(9) 

263.5998326931(2) 

0.733(13) 

220.43(8) 

9.83396979(8) 

9.497573(11) 

54005.80617(6) 

0.053108(3) 

191.500(2) 



17*^50™ 13? 8016(5) 

-37°03'10'.'95(2) 

8.96050624775(3) 

-0.4545(6) 

233.82(3) 

17.3342759(7) 

24.39312(8) 

54003.127812(11) 

0.712431(2) 

131.3547(2) 

0.00548(30) 



17*^50™ 12^1770(4) 

-37°03'22'.'93(2) 

164.62146212170(16) 

-0.520(11) 

234.391(9) 

3.60511446(5) 

2.865858(13) 

54002.7705(11) 

0.004046(9) 

323.07(11) 



0.000822203(16) 



0.00951180(3) 



0.12 



0.30 



0.0518649(5) 

1.97(15) 

1.65 

0.53 



0.00194450(3) 



0.17 



Calculated assuming a pulsar mass of 1.4 Mq. Where a uj was measured, these limits include 99% of the total probability. 



sars, i.e., there are no unmodeled trends in the TOA 
residuals (for two examples, see Fig. [S]). For this reason, 
we believe that the 1-a uncertainty estimates returned by 
TEMPO are essentially accurate, so they are reported 
directly in Tables [3] and 21 In i i4.31 we test the valid- 
ity of this hypothesis for the measurement of the uj of 
NGC 6440B. 

The large EFAC of NGC 6440C indicates the presence 
of unmodeled effects in its TOA residuals. The cause 
of these is unknown at present, and will be investigated 
elsewhere. 

3.1. Positions 

In TableO we calculate the angular offsets of the pulsar 
positions listed in Tables [3] and |4] relative to the centers 



of their respective GCs listed in Table [H These are dis- 
played graphically in the top panels of Figures [4] and [5l 
The projected distances of the pulsars from the center of 
the GC {r±) were calculated using the most recent GC 
distance estimates listed in Table [5] There is nothing 
unusual about the spatial distribution of the pulsars in 
these GCs. As in 47 Tucanae and most other GCs, all 
pulsars are located with in the half-mass radius of the GC 
(|Camilo fc Rasio 20051) : in the present cases at less than 
0f6 (or 4.5 core radii) from the center. This is not an 
observational bias: the S-band telescope beam still has 
> 50% sensitivity 3' from the center of the GC. Further- 
more, the 820-MIIz observations do not reveal any extra 
pulsars within ~ 8' from the center of the GC. The con- 
gregation of pulsars at the center is a real phenomenon. 
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Fig. 3. — Time of arrival residuals as a function of orbital phase for NGC 6440B {top) and NGC 6440D {bottom). 



TABLE 5 

Pulsar offsets from the centers of 

NGC 6440 AND NGC 6441 



Pulsar da"^'^ 


e^" 


e 


j_ 


ri_ 


{') 


{') 


{') 


(Sc) 


(pc) 


NGC 6440 


A ... -0.0017 


-0.0447 


0.04 


0.34 


0.11 


B ... 0.0601 


-0.0307 


0.07 


0.52 


0.16 


C ... -0.3571 


-0.2801 


0.45 


3.49 


1.08 


D ... -0.2472 


0.4930 


0.55 


4.24 


1.32 


E . . . . 0.0240 


0.1269 


0.13 


0.99 


0.31 


F . . . . -0.0858 


-0.0389 


0.09 


0.72 


0.22 


NGC 6441 


A ... 0.1797 


-0.0993 


0.21 


1.87 


0.81 


B ... -0.1432 


-0.2989 


0.33 


3.01 


1.30 


C ... 0.1111 


-0.0096 


0.11 


1.01 


0.44 


D ... 0.0395 


-0.0229 


0.05 


0.41 


0.18 



^ The uncertainty of these parameters is much 
smaller than the uncertainty of the position of the 
center of the GC (they were computed assuming 
that the center is exactly where indicated in Table 
[2J. This precision is, however, necessary to calculate 
the angular distances between the pulsars.'' Note 
that da = da cos 5, where da is the difference of the 
Right Ascension of the pulsar and the center of the 
GC. 

due to mass segregation in the GC. 

3.2. Period derivatives 
In the third column of Table [SI we list the observed 
period derivatives (Pobs)- Four effects can contribute to 
Pobs (|Phinnev 19931) : 



+ -{ac + acc + apM) , (3) 
c 



where Pint is the pulsar's intrinsic period derivative, 
ace is the line-of-sight component of the acceleration 



of the pulsar caused by the gravitational field of the 
GC, aa is the difference of the accelerations of the GC 
and of the Solar System in the gravitational field of the 
Galaxy, projecte d along the Hne of sight, apM = fi^D 
(jShklovskii 1970 ^ is the centrifugal acceleration due to 
the transverse motion of the pulsar /i, and D is the GC's 
distance from the Earth, listed in Table [H The pulsar 
proper motions ^ are not yet directly measurable, and 
there are no optical measurements of the proper motions 
of the two GCs; so, like Pint, apM is not known but is 
always positive. 

In the fourth column of Table [SI we list the observa- 
tional upper limits for the pulsar acceleration, Cp^niax, 
obtained by re-arranging all the known terms of Eq. [31 
to the right side of the equation: 



acc + c 




-aa. (4) 



We calculated ag using a mass model of the Galaxy 
(jKuiiken &: Gilmore 1989f ) , the estimated distance to the 
GCs, and their Galactic coordinates given in Table [2l 
Since we have no estimates of the proper motion of the 
GCs, we assume from now on that iJ,^D = 0. In all GCs 
where this term is measured, it is small relative to the 
intrinsic period derivatives and even smaller compared 
to the accelerations induced by the GC. 

In the fifth column we list the theoretical maximum 
value of ace for each pulsar's line-of-sight, accmax- 
This is calculated using an analytical mass model 
for the centers of GCs described in the Appendix of 
Freire et al. (2005), using the core radius and central 
velocity dispersion listed in Table [2] Both ap^max for 
the individual pulsars and aocmnK for all lines of sight 
close to the center of the GC are displayed in the bottom 
halves of Figures [4| and [5l 

Using Eq. [U we can calculate upper (and, in the cases 
of NGC 6440A and D, lower) limits for Pint assuming the 
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TABLE 6 
Limits for derived parameters of the NGC 6440 and NGC 6441 pulsars 



Pulsar 



P 

(ms) 



-^obs 

(10-18) 



(10- 



QGCmax 

.^ „-2^ 



Pin 



max ^in 
(10-18) 



:iax -Dn 

(10" G) 



(Gyr) 



NGC 6440 



A ... 


288.6027917197(16) 


399.33(10) 


414.8 


26.30 


424.7 


374.0 


353.2 


331.4 


0.012 


0.011 


B ... 


16.760127219257(5) 


-0.32913(16) 


-5.86 


24.72 


1.05 


(0.0) 


4.24 


(0.0) 


{+^) 


0.25 


C ... 


6.226932720487(2) 


-0.05984(16) 


-2.85 


6.15 


0.07 


(0.0) 


0.66 


(0.0) 


{+^) 


1.45 


D ... 


13.495820403909(2) 


0.58678(16) 


13.06 


4.77 


0.80 


0.37 


3.32 


2.26 


0.57 


0.27 


E ... 


16.264003411125(5) 


0.3124(4) 


5.79 


19.62 


1.38 


(0.0) 


4.78 


(0.0) 


(+oo) 


0.19 


F ... 


3.793629114949(3) 


-0.01055(18) 


-0.81 


22.53 


0.27 


(0.0) 


1.03 


(0.0) 


(+^) 


0.22 



NGC 6441 


A ... 111.6008373133(4) 
B ... 6.074542086503(6) 
C ... 26.56867837497(4) 
D ... 5.139936774902(4) 


5.661(8) 
0.0192(4) 
-0.996(3) 
0.4928(2) 


15.67 

1.41 

-10.77 

29.20 


19.88 
11.88 
31.30 
41.46 


13.2 
0.27 
1.82 

1.21 


(0.0) 
(0.0) 
(0.0) 
(0.0) 


38.8 
1.29 
7.01 
2.52 


(0.0) 
(0.0) 
(0.0) 
(0.0) 


( + 00) 

(+oo) 
(+oo) 


0.13 
0.36 
0.23 
0.067 



extreme possible accelerations iaccmax- These limits 
are displayed in the sixth and seventh columns of Table 
[SI From these Pint limits, we can calculate upper and 
lower limits for the magnetic field at the surface (esti- 
mated using Bo = 3.2 x lQ^^{PPi^tYl'^ G) and lower 
and upper limits for the characteristic age of these pul- 
sars (estimated using t^ — P/2Pi^^t), these are displayed 
in the final columns of Table [HI No pulsar in either GC 
has Upjnax < — OGCmax, this means that the present GC 
mass models can predict line-of-sight accelerations large 
enough to explain the observed negative period deriva- 
tives (and plausibly maintain that the pulsars are not 
spinning up). 

3.2.1. NGC 6440A 

As previously found by Lyne, Manchester & D'Amico 
(1996), PSR B1745-20A has a characteristic age ~ 10^ 
times smaller than the age of the GC and Bq ^ 3 x 
lO^^G, a value that is more typical of what one finds in 
the general Galactic population. For these two reasons, 
there have been some doubts about the association of this 
pulsar with NGC 6440. The discovery of five new MSPs 
in NGC 6440 at DMs similar to that of PSR B1745-20A 
confirms its association with NGC 6440. Furthermore, 
PSR B1745-20A is closer to the center of the GC {6^ = 
0f04) than any of the new MSPs. 

This association of a "young" pulsar with a GC is 
not unique to NGC 6440. The first known examples 
were PSR B1718-19, in NGC 6342 (iLvne et al. 19931) 
and PSR B1820-30B in NGC 6624 (iBiggs et al. 19941) . 
The possible origins of these objects, and in particular 
NGC 6440A, are discussed in detail in Lyne, Manchester 
& D'Amico (1996). Either neutron star formation contin- 
ues to happen in GCs — through electron capture super- 
novae (ECS), which can form when an accreting WD star 
reaches the Chandrasekhar limit or when a WD binary 
coalesces — or some form of mild recycling happened, 
in which the pulsar was not spun up to very short spin 
periods and not much of its magnetic field was buried. 
The ECS scenario is particularly appealing as this type of 
supernova also produces small kicks, which nicely solves 
the problem of neutron star retention in GCs (Ivanova 
et al. 2007). 

3.2.2. NGC 6440D 

Unlike in the case of NGC 6440A, the upper limit for Tc 
and the lower limit for B does not imply that this pulsar 



is younger or has a stronger magnetic field than the other 
MSPs. It is possible that the other "slow MSPs" (with 
13 < P < 27 ms) have similar ages and magnetic fields 
as those of NGC 6440D; their intrinsic Ps might simply 
being masked by the GC accelerations. The apparently 
special limits for NGC 6440D might have more to do with 
its large distance from the center of the GC: the range of 
possible GC accelerations aocmax is the smallest for all 
pulsar positions in NGC 6440 and NGC 6441; this leads 
to tighter limits on Pint- If NGC 6440D is a typical slow 
MSP, then the magnetic fields of these objects are one 
order o f magnitude large r than those observed in 47 Tu- 
canae ()Freire et al. 200 ll) ; this would cause the observed 
difference in spin periods between GCs discussed in ij2.2l 

4. BINARY SYSTEMS 

Before this survey, only one binary pulsar was kno wn in 
these two GCs, NGC 6441 A (jPossenti et al. 20061) . We 
have discovered four new binary systems, three in NGC 
6440 (B, D and F) and one in NGC 6441 (B). These bina- 
ries make up half of the total observed pulsar population 
in these GCs. This ratio is similar to that of Terzan 5, 
but significantly lower than what we find in 47 Tucanae 
and M62, in the latter GC all the 6 known pulsars are in 
binary systems. 

4.1. NGC 6440F and NGC 6441B 

These two binaries, with orbital periods of 9.8 and 3.6 d 
and minimum companion masses of 0.30 and 0.17 Mq 
respectively (assuming pulsar masses of 1.4 Mq), might 
be the only "normal" MSP- WD binaries among the new 
pulsars. Their spin periods (3.79 and 6.07 ms) are at 
the lower end of the spin period distribution found in 
NGC 6440 and NGC 6441, but they arc typical of what 
is found among similar MSP- WD binaries in other GCs 
and in the Galactic disk. Their eccentricities (0.0531 
and 0.00404) are small compared to the highly eccentric 
NGC 6440B and NGC 6441A, but they are nevertheless 
much larger than what is found in binary systems with 
similar orbital periods in the Galactic disk. This is a 
common occurrence in GCs, and is probably due to per- 
turba tions caused by close flybys of other stars in these 
GCs (jHeggie fc R.asio 19961 ). 

4.2. NGC 6440D - Eclipsing binary 

The binary system NGC 6440D has an orbit with non- 
measurable eccentricity and a period of 6.9 hr. It is an 
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Fig. 4. — Top: Positions of the six known pulsars in NGC 6440, 
plotted as west-east {6a) and south-north (8g) offsets from the 
center of the GC (see also Table [Sjl. The angular core radius {9c 
= 0fl3) is represented by a circle. Bottom: Upper limits for the 
accelerations of the MSPs {triangles) . These are compared with 
the maximum acceleration that can be caused by the GC, as a 
function of 6± {solid curves). The angular core radius is indicated 
by the vertical dashed line. Pulsar NGC 6440A, and to a much 
lesser extent NGC 6440D are significantly above the maximum 
GC acceleration, indicating that its large positive P is caused by 
its intrinsic slowdown. 



eclipsing system with eclipses that usually last for about 
10% of the orbital cycle. Its minimum companion mass 
is 0.12 Mq, calculated assuming a pulsar mass of 1.4 Mq. 
Freire (2005) called attention to the fact that eclips- 
ing binary pulsars have a bi-modal distribution of mass 
functions: those with / > 1O~^M0 were designated 
"Eclipsing Low-Mass Binary Pulsars" (ELMBPs), while 
those with / < 3xl0~^ Mq were designated as "Very 
Low-Mass Binary Pulsars" (VLMBPs). The VLMBPs 
are more commonly known as "Black Widow" pulsars. 
NGC 6440D is clearly a ELMBP. Unlike "Black Wid- 
ows", these objects only occur in GCs, suggesting that 
they form through exchange interactions when a radio 
pulsar acquires a MS companion. 

Sometimes NGC 6440D is detectable during superior 







arcmin 



Fig. 5. — Top: Same as for the top of Fig.|4]for the four MSPs in 
NGC 6441. The angular core radius {6c = O.'ll) is represented by 
a circle. Bottom: Upper limits on the accelerations of the MSPs in 
NGC 6441 (triangles) and theoretical maximum acceleration as a 
function of 6± (solid curves). 



conjunction. This indicates that the orbital inclination 
is certainly less than 90° , otherwise the companion itself 
would cause an eclipse. It also implies that the concen- 
tration of the material producing the eclipse is highly 
variable with time. Nevertheless, the cloud of material 
seems to generally lead the companion in orbital phase 
(see Fig. El, i.e., ingress occurs at a much larger distance 
from superior conjunction than egress. Before, after, and 
sometimes through the eclipse we detect '--^O.Sms delays 
in the times of arrival (see Fig. [3]). These are proba- 
bly due to an increase in the electron column density 
near the companion of up to ~0.5cm^'^pc (or about 
1.5 X lO^^cm^^). This is of the same order of magni- 
tude as wh at has been found in similar binary systems, 
e.g. M30A (jRansom et al. 2004D . 

If the system is nearly edge on, with pulsar and com- 
panion masses of 1.4 and 0.124 Mq, the separation be- 
tween the stellar components is 1.47x10^^ cm {2.1 Rq). 
On some occasions the eclipse extends up to ^ 30° from 
superior conjunction (normally before the companion). 
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which implies the detection of material at ^ 1.0 Rq from 
the center of the companion. The size of the Roche lobe 
of the companion is 0.42 Rq, therefore, it is clear that we 
are detecting the presence of material that is not bound 
to the companion at ingress. We almost never detect 
unbound material at egress. 

This is different fro m what one observ es in ELMBPs 
such as 47 Tuc W dC amilo et al. 2000[ ), Terzan 5 P 
and ad (| Ransom et al. 2005; Hcss els et al. 2006I) . 
PSR J1740-5340 in NGC 6397 (|D'Amico et al. 2001h 
or M28 H (Begin 2006; Begin et al. in preparation), 
which display long duration irregular eclipses likely 
related to excessive unbound ionized material from 
the companion over much of the orbit. Three of these 
ELMBP companions have been identified at optical 
wavelengths (47 Tuc W, Edmonds et al. 2002, ; 
PSR B1718-19 in NGC 6349, Kerkwijk et al. 2000; and 
PSR J1740-5340 in NGC 6397, Ferraro et al. 2001 ), 
and in all cases the observed companion is a low-mass, 
non-degenerate star. None of the EL MBPs with shorter 
eclipses such as M30 A, Terzan 5 A (jLvne et al. 19901) 
and M62 B (jPossenti et al. 2003f ) have been optically 
detected. The situation is similar at X-rays; only the 
systems with large eclipses have been detected (e.g., 
47 Tuc W, see Bogdanov, Grindlay & van den Berg 
2006). In these systems the hard X-ray emission is gen- 
erated by the collision of MSP and stellar winds. None of 
the ELMBPs with shorter eclipses have been definitively 
detected at X-ray wavelengths. If NGC 6440D is similar 
to these ELMBPs, it seems less likely that it will be 
detected at either optical or X-ray wavelengths. 

When an exchange encounter forms a ELMBP, the 
ejection of the previous companion to the pulsar imparts 
a kick to the binary, making it fly away from the central 
regions of the GC in an eccentric orbit. It is probably 
for that reason that two ELMBPs (PSR B1718-19 and 
PSR J1740— 5340) are found at large dist ances from the 
centers of their parent GCs (jFreire 2005D . NGC 6440D 
follows this trend; of all the pulsars in NGC 6440 and 
NGC 6441, it is farthest from the center of its GC as 
projected on the plane of the sky. 

4.3. NGC 644OB: eccentric binary with a super-massive 
neutron star? 

NGC 6440B and NGC 6441A {SM are part of a di- 
verse class of binary pulsars in GCs with eccentric or- 
bits that includes m ostly recent discoveries like M30B 
(jRansom et al. 20041 ). NGC 1851A (Freire et al. 2004; 
Freire et al. 2007), six systems in Terzan 5 (I, J, Q, U, X 
and Z; Ransom et al. 2005; Stairs et al., in preparation) 
and M28 C and D (Begin 2006; Be gin et al. in prepara - 
tion). With the exception of M15C (jjacobv et al. 2006f ). 
which is very similar to some DNSs seen in the Galaxy, 
they have no Galactic analogues. All of these systems 
were likely formed through stellar exchange encounters. 

The orbital eccentricity of NGC 6440B has allowed a 
highly significant measurement of the rate of advance 
of periastron; ui = 0.00391(18)°yr^^. As mentioned 
above, this is a l-cr uncertainty derived by TEMPO. To 
check if this uncertainty value is realistic, we estimated 
to and its associated uncertainty using two other meth- 
ods. First, all the timing parameters of NGC 6440B 
were estimated using a Monte-Carlo Bootstrap algo- 
rithm (jEfron fc Tibshirani 199llPress et al. 19921) . with 



a total of 1024 fake TOA datasets that are consis- 
tent with the real dataset. From this we obtain u; — 
0.00391 (16)°yr~^. In the second estimate, we kept uj 
fixed and fitted all the remaining timing parameters, 
recording the resulting x^- Doing this for a range of val- 
ues of oj, we can estimate the l-cr uncertainty as the half- 
width of the region where [x^ (w) — Y^(a;min) < IL bein g 
Winin the value that minimizes x^ (jSplaver et al. 200^ . 
The result of this estimate is w = 0.00391(18)°yr-i. 
Both these estimates indicate that, at least for the u! of 
NGC 6440B, the 1-ct TEMPO uncertainty is rehable. 

The orbital coverage of the timing data is excellent (see 
Fig. [3|); the orbit was well sampled at the end of 2005, 
at the start of 2007 and again during October 2007. For 
this reason, tb is not strongly correlated to any timing 
parameter: its largest correlation is with Pt, (96%), which 
is normal for this particular parameter. 

4.3.1. Mass estimate from periastron advance 

Assuming that uj is fully relativistic (an important as- 
sumption, as we will see) implies a total system mass of 
2.92 ± 0.20 Mq. This binary system is located, at least in 
projection, close to the center of the GC, where the mas- 
sive objects should predominantly occur. In Figure[6l we 
present the mass constraints on the components graph- 
ically. We assume that the probability density function 
(pdf) for UJ, p{uj), is a Gaussian, with half- width similar 
to the 1 cr uncertainty given by TEMPO. We also as- 
sume a constant probability density for cosi, where i is 
the orbital inclination (90° for an edge-on orbit). From 
these constraints, we calculate a two-dimensional pdf for 
the mass of the pulsar and the mass of the companion. 
This is then projected in both dimensions, resulting in 
final pdfs for the mass of the pulsar and the mass of the 
companion. 

The total mass of this system suggests that it is a 
double neutron star (DNS) binary. However, the prob- 
ability that the pulsar lies within the range of neutron 
star masses that have been precisely measured to date — 
from ~ 1.20 Mq for t he companion of PSR J1756-2251 
(jFaulkner et al. 20051) to 1.44 Mq for PSR B1913+16 
(jWeisberg fc Tavlor 2003^ , indicated in Figures [5] and [7] 
by grey bars — is only 0.10%. This eventuality would 
require very low orbital inclinations, between ~ 4° and 
^5°. In this inclination range, it is also possible that 
the pulsar has a blue straggler companion, a possibility 
that can be investigated using optical/IR imaging. 

Before this study, the most massive neutron star known 
was PSR B1516-h02B in the GC M5 (1.96+[j;?^M0; Freire 
et al. 2007). If NGC 6440B had a similar mass, then the 
companion could be a ^ 0.9 Mq white dwarf, or a blue 
straggler of similar mass. However, we note that there is 
only a 0.97% probability that the mass of the pulsar is 
below 2.0 Mq. 

The median pulsar mass is 2.74 Mq, with the lower 
and upper la limits at 2.52 and 2.95 Mq, and the 2cr 
limits at 2.23 and 3.15 Mq. The pulsar has a 99% prob- 
ability of being less massive than 3.24 Mq. For the com- 
panion, the median of the distribution is at 0.142 Mq, 
with the the lower and upper la limits at 0.124 and 
0.228 Mq, and the 2 a limits at 0.113 and 0.571 Mq, im- 
plying that it could be a low-mass WD or an un-evolvcd 
MS star. If this high pulsar mass is confirmed, it would 
be by far the largest neutron star mass ever measured. 
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Fig. 6. — Graphical display of the constraints on the masses of NGC 6440B and its companion. In the main display, the hatched region 
is excluded by knowledge of the mass function and by sini < 1. The slanting straight lines correspond to a total system mass that causes a 
general-relativistic ili equal or within 1 tr {dashed lines) of the measured value. The five solid curves indicate constant orbital inclinations. 
The gray bars indicates the range of precisely measured neutron star masses. We also display the probability density function for the mass 
of the pulsar (top) and the mass of the companion (right), and mark the respective medians with vertical (horizontal) lines. 



That would have profound consequences for the study of 
the equation of state of dense matter, since almost no 
models predict neutron sta rs more massive than 2.5 Mq 
(jLattimer fc Prakash 2007D . 

If NGC 6440B really has a mass of - 2.7 A/©, it 
would imply that the end products of the coalescence 
of DNS systems might themselves be stable as super- 
massive neutron stars. The total masses of the known 
DNS system s range from 2.57 M^^ (iFaulkner et al. 2005( 1 
to 2.83 Mq (jWeisberg fc Tavlor 20031 ) . Such coalescence 
products might be observable after c oalescence, eithe r 
through gravitational wave emission (jAndersson 2003f ). 
cooling through neutrino emission (direct URCA process. 
Page & Applegate 1992) or even high-energy blackbody 
emission. 

It is possible that NGC 6440B itself formed this 
way. M15C, the o nly DNS known in the GC system 
(| Jacobv et al. 20061 ) . is expected to coalesce in about 
3 X 10*^ years (^ 3% of the age of M15). If the prod- 
ucts of DNS coalescence are as likely to be detected as 
radio pulsars as their progenitors, then we should detect 
~ 30 of them in the whole GC system. It is possible 
that some isolated MSPs in GCs and the Galactic disk 



formed this way. The massive isolated MSPs in GCs can 
later acquire a companion through exchange encounters, 
forming a system like NGC 6440B where we can measure 
their large masses. 

4.3.2. Is the periastron advance relativistic? 

The previous discussion rests on the assumption that 
the observed to is fully relativistic. If there was an extra 
contribution to io from tidal or rotational deformation of 
the companion, then no reliable estimates can be made 
of the total mass of the system. 

Lai, Bildsten fc Kaspi (1995), in their analysis of the 
binary pulsar PSR J0045— 7319, concluded that the only 
likely contribution to ld in that system is from rotational 
deformation of the companion. NGC 6440B has a shorter 
orbital period than PSR J0045-7319 (20.5 versus 51 
days), but it also has a much smaller companion mass 
(- 0.12 Mq at i = 90° to ~ 1.5 Mq at i ~ 4 - 5°) 
versus ^ 9 Mp, f or the companion of PSR J0045-7319 
(|Bell et al. 19951 ). Main-sequence stars of such masses 
have radii o f R ^ 0.18 — 1.5 and ~ 4.5 i?Q respectively 
(jLang 1991[ ). This implies that even if the companion 
of NGC 6440B were extended, the contribution to uj 
from its tidal deformation (the latter proportional to 
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(R/a)^, where a is the separation between components) 
would always be smaller for NGC 6440B than it is for 
PSR J0045-7319. 

The only possible exception to this is if the compan- 
ion of NGC 6440B were a ~ 0.9 M© star undergoing a 
pre-giant or giant phase. Given the mass function of 
the system, that would require a very low orbital incli- 
nation. A giant star of that mass would have a radius 
of up to 1 a.u. (jLang 1991[ ). and would not fit in the 
space between the pulsar and companion at any orbital 
phase. If it were in a "pre-giant" phase, its atmosphere 
would be very extended and have significant mass loss, 
leading to variations in the DM with orbital phase. The 
resulting tides should not only have circularized the or- 
bit, but they should also have directly observable effects 
on the timing: unmodeled effects in the rotation of the 
pulsar and random variations of the orbit's period and 
apparent size (Nice, Arzoumanian & Thorsett 2000). For 
PSR J0045 — 7319, there is strong orb ital decay due to 
tidal effects (jKumar fc Quataert 1998f ). 

For NGC 6440B, we obtain A = (-0.0 ± 1.4) x 10^^ 
and X = (-0.21 ± 0.14) x 10"^^. The pulsar times very 
well, with an EFAC similar to those of the isolated MSPs 
in the GC (see Table |3]), no unmodeled systematic trends 
are visible in the TOA residuals (Fig. [3]). Fitting for the 
halves of the orbit "in front" and "behind" the compan- 
ion, we see no DM variations larger than our detection 
limit of 0.07 cm~^pc. The possibility of an unusually 
extended companion to NGC 6440B and of a significant 
tidal contribution to cj can be safely excluded. 

The contribution from rotational deformation, hence- 
forth designated as Wrot, could be significant if the 
companion (degenerate or not) were spinning rapidly. 
The rotational deformation will also induce a change 
in the projected semi-major axis of the orbit x. 
Splaver et al. (2002), based on the results of Lai, Bildsten 
& Kaspi (1995) and Wex (1998), relate Wj-ot to x: 



X 1 

w,ot = - tani^— . 

X \ sm t^ cos w sm $0 



I sin^ ( 



cot $0 



(5) 



where i is the orbital inclination of the binary, 6 is the 
angle between the angular momentum vector of the sec- 
ondary and the angular momentum vector of the orbit; 
and $0 is the longitude of the ascending node in a refer- 
ence frame defined by the total angular momentum vec- 
tor (see Fig. 9 of Wex 1998). 

Thus, our observed upper hmit of \x/x\ < 8 x 10"^*^ 
implies |wrot| < 1-5 x 10~'*°yr~^ times a geometric fac- 
tor. In 80% of cases this geometric factor will be smaller 
than 10. Therefore Wrot will in most cases be of the same 
order of magnitude as the present measurement uncer- 
tainty for uj. Certain special combinations of i, 0, and 
$0 will make the geometric factor large, but, again, this 
would only occur if the companion were rotating fast. 

To summarize, the contribution to the uj of this system 
from the tidal deformation of the companion is negligible. 
Furthermore, observational limits on x indicate that uj^ot 
is also likely to be negligible and the observed uj due to 
the effects of general relativity. 

4.4. NGC 64^4^1 A: eccentric binary with a massive 
companion 



NGC 6441 A was the only puls ar previously known 
in NGC 6441 (jPossenti et al. 20061 ). ahhough its phase- 
connected timing solution has never been published. Its 
orbital period is 17.3 days, and its orbit is highly eccen- 
tric (e — 0.71), allowing a measurement of the advance 
of periastron for this system, tj = 0.00548(30)°yr~^. As- 
suming, again, that this effect is fully due to general rel- 
ativity, we obtain a total mass of 1.97 ± 0.15 Mq, con- 
sistent with the estimate made by Possenti et al. (2006): 
2.15 ± 0.06 Mq. Using the methods applied in the case 
of NGC 6440B, we calculated pdfs for the mass of this 
pulsar and its companion for our measurement of the to- 
tal mass of the binary; these are displayed graphically in 
Figure [71 There is a 99% probability that the compan- 
ion is more massive than 0.53 Af© and that the pulsar 
is less massive than 1.65 M©. The medians for the pul- 
sar and companion masses are 1.26 and 0.67 Mq respec- 
tively, and there is a 45.6% probability that the pulsar 
mass is within the 1.20-1.44 Mq mass range. The na- 
ture of the companion is unknown, but no observational 
hints of tidal effects (i, Pb, DM variations) are present, 
suggesting it is a compact object. 

5. CONCLUSION AND PROSPECTS 

We have discovered eight new pulsars in the GCs 
NGC 6440 and NGC 6441. Their pseudo-luminosities 
indicate that there may be as many pulsars in these GCs 
as in Terzan 5, but less pulsars are observed in the for- 
mer because of their larger distances. The pulsar popu- 
lation in these GCs seems to have distinctly lower spin 
frequencies. To some extent, this is due to the presence 
of apparently young objects like NGC 6440A. At least in 
the case of NGC 6440A and D, and possibly in the case 
of the other slow MSPs, the larger spin periods could be 
related with many of them having higher magnetic fields 
than the MSPs in other GCs. 

Four of the new pulsars are in binary systems. 
NGC 6440D is an eclipsing binary, with a companion star 
with a mass of about 0.12 M© or larger. NGC 6440B is a 
16.7-ms pulsar in a 20.5-day orbit with an orbital eccen- 
tricity of 0.57. A measurement of the rate of advance of 
periastron for this system suggests that this binary hosts 
the most massive neutron star to date: 2.74 ± 0.21 Af©, 
with only a 1% probability that the inclination is low 
enough that pulsar mass is below 2 Af © . This result de- 
pends on the (likely) possibility that the observed peri- 
astron advance is fully relativistic. Such a large mass 
would introduce the strongest constraints to date on the 
equation of state; it also suggests that the products of 
the coalescence of double neutron stars might be stable 
neutron stars. The rate of advance of periastron was also 
measured for NGC 6441A, indicating a significantly less 
massive pulsar {Mp < 1.65 M©). 

The pulsars in NGC 6440 will be monitored carefully 
over the next few years. In the case of NGC 6440C, the 
new data will allow a detailed investigation of the origin 
of its timing irregularities. In the case of NGC 6440B, 
the continued monitoring will lead to an improvement in 
the precision of Co. Simulations indicate that the Ein- 
stein delay (7) may become measurable with about 15 
years of data; this would provide unambiguous mass mea- 
surements and determine the nature of the companion. 
Deep optical studies of NGC 6440 might be extremely 
valuable in determining the nature of the companion of 
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Fig. 7. — Display as in Fig. (6] this time for the NGC 6441A binary system. 



NGC 6440B. 
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